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SUMMARY 

The equations governing the three-dimensional compressible 
laminar boundary layer with variable viscosity and thermal con- 
ductivity are shown to be simplified by: 

1. The introduction of a two -component vector potential 

2. The use of a transformation to change the equations into 

nearly incompressible form 

3. The use of a further transformation that changes the equa- 

tions into nearly Cartesian form when a coordinate system 
appropriate to axial or conical symmetry is used 

Flow over flat plates with arbitrary leading-edge contours is 
discussed and it is deduced that, under certain circumstances, the 
boundary -layer equations are inapplicable in certain restricted 
regions of the boundary layer* 

Problems involving conical potential flow are discussed, and 
it is shown that use of the Blasius similarity variable permits 
reduction of the number of independent variables. 


INTRODUCTION 

The techniques for predicting the behavior of the two-dimensional 
laminar flow of a viscous fluid are rather well developed. Because 
of the interest in aircraft applications, most of the effort has been 
applied to problems of the boundary layer, which is associated with 
flight at substantial speeds through a medium of low kinematic 
viscosity. An investigation of methods for treating the three- 
dimensional compressible laminar boundary layer was conducted at the 
NACA Lewis laboratory and is presented herein. 

Certain features of the present body of theory for plane cases 
furnish guidance: 
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1. Lagrange's stream function nay be introduced in order iden- 
tically to satisfy the equation of continuity, thus permitting two 
of the dependent variables (velocity components) to be expressed in 
terms of a single f unction. 

2. In the absence of a pressure gradient, Blasius (reference 1, 
paragraph 53) shows that, for the incompressible case, the two inde- 
pendent variables (space coordinates) may be combined into a single 
coordinate of similarity, thus reducing the problem to that of solving 
an ordinary differential equation for the stream function. 

3. The Barman -Pohlhausen integral method (reference 1, para- 
graph 60) provides a valuable engineering approach, the simplicity 
of which is obtained at the expense of restriction to a particular 
family of velocity profiles. 

4. Howarth (reference 2) shows that a certain transformation 
of coordinates results in momentum equations of nearly incompres- 
sible form, provided that a linear dependence of viscosity on tem- 
perature is assumed. 

5. With regard to the energy equation, if the Prandtl number 
is 1 and the wall is an insulating surface, the enthalpy remains 
constant through the boundary layer. Furthermore, if no pressure 
gradient exists, a term in the solution for temperature depending 
linearly on the velocity profile may be introduced in order to take 
into account heat transfer from a wall at constant temperature 
(Crocco, reference 3, and von Barman, reference 4). 

In all the theories previously mentioned and in the present 
report, use is made of the usual order -of -magnitude assumptions pro- 
ceeding from the concept of a thin boundary layer beyond which exists 
potential flow undisturbed by the presence of the boundary layer. 

(See, for example, reference 1, paragraphs 44 and 45.) 

The three-dimensional problem is, of course, greatly complicated 
by the number of dependent and independent variables appearing in the 
equations of motion. When such difficulty is encountered, an inte- 
gral approach would appear to provide adequate simplification of the 
problem. Prandtl (reference 5) proposes a method (for incompressible 
flow) utilizing two parameters; namely, a boundary -layer "thickness" 
and the local angle of divergence between the outer streamlines and 
the limiting streamline at the wall. This method is not the only 
possible way of formulating an integral approach; in fact, the very 
complexity of the three-dimensional problem tends to make it possible 
to devise a number of such formulations. For example, the two 
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independent parameters can "be defined as the two boundary-layer thick- 
nesses S u and corresponding to the two velocity components u 

and w in planes parallel to the -wall. This formulation would 
presumably be useful if it were expected that one component would 
approach its stream value much more rapidly than the other. 

The relative merits of Prandtl' s formulation and some other, such 
as the one just described, are difficult to assess; nor is it a 
simple matter to make a judicious selection of profile functions to 
be used in a three-dimensional problem. Prandtl (reference 5) points 
out that the setting-up of an integral method requires a background 
of experience gained from solutions to the complete equations of 
motion, lor this reason, the present report is concerned with means 
of attacking the differential boundary -layer equation and will not 
discuss integral methods further. 

The equations of motion have been solved in several three- 
dimensional cases: Prandtl (reference 5), Sears (reference 6), and 

R. T. Jones (reference 7) have contributed to the incompressible solu- 
tion for a yawed infinite cylinder. The most striking result in this 
problem is that the velocity components in a plane normal to the 
cylinder may be obtained by two-dimensional theory and may then be 
used to determine the axial component. Problems involving spherical 
or axial symmetry have received considerable attention. Mangier 
(reference 8) provides an analysis of the axially symmetric incompres- 
sible flow over bodies of revolution, wherein various changes of vari- 
able permit the direct use of plane boundary- layer theory. In unavaila- 
ble work. Mangier has extended this analysis to permit consideration of 
compressibility. 

The solutions mentioned in the preceding paragraph are obtained by 
discovering methods by which plane-flow results may be directly adapted. 
Many problems in the mechanics of the laminar boundary layer remain in 
which this simplification is impossible and the analysis to follow will 
investigate methods of solving this type of problem. . 

SYMBOLS 

The following symbols are used in this report: 

C constant appearing in viscosity -temperature relation 
Cp specific heat at constant pressure 
c v specific heat at constant volume 
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E internal energy, c v T 

F component of vector potential in implicit coordinates 

f component of vector potential in conical coordinates 

G component of vector potential in implicit coordinates 

g component of vector potential in conical coordinates 

1 o P 

3 total enthalpy in boundary layer, cpT + JjIu +w ) 

k arbitrary constant 

Pr Prandtl number , Cpp/tt 


P 

R 

r(x) 

S=B 


T 

t 


u,v,w 


x,y,z 


I,Y,Z 


5 

9 

K 


pressure 
gas constant 

distance used in definition of implicit coordinates 
dimensionless coordinates 
absolute temperature 
time 

velocity components 
coordinates 


equal to x. 



fjr 


_p_ 

p« 


dy, and z, respectively 


>0 


boundary -layer thickness 
enthalpy, c^T 

coefficient of thermal conductivity 
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_1 

2 

\ Blasius similarity variable, t] i 

U coefficient of viscosity 


I »T| , | equal to k 


X 


r ax, krY, and S, respectively 


P 

a 


density 

similarity variable for flat plate. 


1 

2 


Y(X-X 0 ) 


cp,cp 


components of vector potential 


Subscripts j 

1 denotes conditions in flow at outer edge of boundary layer 

» denotes evaluation at some reference condition 

0 denotes conditions at leading edge (See figs. 1 and 2.) 

Subscript notation for partial differentiation is used where convenient. 
Primes denote ordinary differentiation. 


THEORY 

Equations of Laminar Compressible Boundary Layer 

In Cartesian tensor notation, the equations governing the motion 
of a compressible viscous gas may be written as 
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(sr + u ^) ' ‘ + p ( 


r duj_ du 

5^ + 



^ + 4 (PU3> ' 0 


0 dE , dE 

p 5t + 




dui 


>(D 


Sxj 


a 

SxJ 


tij) 


p « PET 


- Equations of motion in Cartesian coordinates* - On the assumption 
of " "boundary -layer flow” an obvious extension to three dimensions of 
the argument presented in reference 1 (paragraph 44) permits equa- 
tions (1) to he reduced to the following form , where a Cartesian 
system of coordinates (x,y,z) and velocity (u,v,w) is used; y and 
v are taken normal to a plane surface hounding the flow: 


p fei * v&i + v^i * w^ - - h + fiffl 

\ot ox oy ozj ox oy \ oyy 


(&) 


|e *o 


(2h) 




_ + a 

5r + 3y 



(2c) 


rr •+ -£-(pu) + -IUp^) + -r-(P w ) “ 0 

at ox oy dz 


(2d) 
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p = PRT (2f ) 


The argument contained in reference 1 (paragraph 45) indicates 
that for plane flow, these Cartesian equations apply in general 
orthogonal coordinates, where y is taken normal to the surface, and 
thus apply to flow over curved surfaces, provided that the surface 
curvature is moderate and with the exception that, for curved surfaces^ 
dp /ay is, in general, of order 1. Of course, even through dp/dy is 
of order 1, the change in pressure across a thin boundary layer is 
small and may he neglected. This argument applies equally well to 
three-dimensional flows. 

Equations of motion in implicit coordinates. - The boundary -layer 
equations are now written in an orthogonal coordinate system wherein: 

1. The body surface is defined by y » 0. 

2. A point is defined by the distances x, y, and r(x) s 
where the distance r(x) s depends implicitly on the distance x, 
and where r(x) has the dimensions of length. 

This coordinate system would be useful in the analysis of flow 
about bodies for which a coordinate x can be defined such that body 
cross sections are similar for various values of x. The quantity 
r(x) then gives the variation of scale of these cross sections. For 
a body of revolution, x may be measured along generators ; the cross 
sections are circular; and r(x) my be taken as the radius of a 
cross section. For a conical body, x my be measured along rays 
from the apex, and r(x) is a linear function of x, giving the scale 
change of the (in general) noncircular cross sections. The velocity 
components are taken to h© u, v, and w in the directions x, y, 
and s, respectively. The equations of motion are : 


p + v|H + w± ^ - El w 2 ) * - & + 

\ot ox oy r os r / dx oy V oy) 


(3a) 
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dp = order of 1 

dy 

(3b) 

p + Y&L + wi + El - - i & + 

\dt ox ay r os r J r ds ay ^ ojj 

| (3c) 

d b b 1 d 

^(pr) + ^(pru) + ^(pry) + - ^(prw) = 0 

(3d) 

_ /d E dE dE 1 dE \ id/ % dv , 1 dw 

p V5t + + v 5y + T w 5sJ = " p [r SI (ru) + SJ + ? 

+ 

/ du\ 2 , ( dw\ 2 . d ( dT^ 

p [_\3y/ W J + V k 5F; 

(3e) 


p - PET (3f ) 


Simplification of energy equation. - Equations (2d) and (2f ) may 
Tse combined with equation (2e) to yield 



'where e s CpT* The specific heats c v and Cp, and the Prandtl 
number Pr = °pp / k are considered constant. Furthermore, if the 
Prandtl number is assumed equal to 1, the sum of equation (2a) mul- 
tiplied hy u, equation (2c) multiplied by w, and equation (4) 
yields, for steady flow, 

_ ( dH dH 8 h\ d ( dll\ 

P + V 5F + V ZZJ = Zy V^J 


( 5 ) 
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■where 


H z CpT + |(u 2 +v 2 ) 


( 6 ) 


A special solution of equation (5) is, in accordance with defi- 
nition (6), 


H * CpT + i(u 2 +w 2 ) 


constant 


(7) 


This solution is presented in references 3 and 4 for the plane case. 
Evaluation of the derivative of equation (7) at the wall, where 
u ■ w ■ 0, shows that equation (7) is the solution of the energy 
equation (5) for the case of zero heat transfer through the wall 
y = 0. The result of equation (7) is, of course, independent of 
the coordinate system used and hence is a special solution of 
equation (3e). 


when 


In references 3 and 4 
Pr = 1 and 


1 = °’ 


, it was observed that for plane flow 

j * 

a solution to the energy equation is 


> 


H = A + Bu (8) 

where A and B are constants to be determined from the boundary 
conditions on the temperature profile. In three dimensions, it is 

clear that if ^2 (but not necessarily ^2) vanishes, then the steady- 

state form of equation (2a) for u is identical to equation (5) for H* 
Thus, equation (8) is a solution of equation (5). Inasmuch as the 
velocity components vanish at the wall, equation (8) must correspond 
to cases involving heat transfer from a wall at constant temperature. 
Because H must be constant in the outer potential flow, solutibn (8) 
is further restricted to cases wherein uq , (the value of u at 

outer edge of boundary layer) is constant. An example of a situation 
wherein solution (8) applies is furnished by the yawed infinite cyl- 
inder, where u is in the direction of the cylinder axis. The anal- 
ogous result for incompressible flow is mentioned in reference 6. 


. Vector Potential 

Definitions . - In order to reduce the number of dependent vari- 
ables appearing in the boundary -layer equations, it is customary in 
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problems of steady plane flow to Introduce the stream f unction, 
which identically satisfies the equation of continuity (equation (2d))j 
that is. 


Pu 


at 

H Sy 


Pv * 


at 

5x 


a 

> O) 

J 


It would be desirable to secure a similar advantage in the analysis 
of three-dimensional steady flows, which may be done by writing 


PU 



"\ 


Pv 


ay _ ap 

5x 5z 


> ( 10 ) 


pw = 


acp 

¥ 


J 


Equations (10) represent one of several ways in which two functions 
can be defined in order to satisfy equation (2d) identically for 
steady motion. The particular arrangement of equations (10) is chosen 
to provide symmetry of y and x against 9 and z. Of course, for 
plane flow, equations (10) reduce to equations (9). 

For every steady plane flow satisfying the equation of continuity, 
a stream function exists, according to definition (9). Proof of this 
theorem is provided by Lamb (reference 9, paragraph 59). Existence 
of the functions cp and y defined by equations (10 ) for every 
three-dimensional flow should, however, be proved. 

A well-known theorem of vector analysis states that any contin- 
uously differentiable vector having a vanishing divergence may be 
expressed as the. curl of a vector potential. By the equation of con- 
tinuity, the mass-flow vector p<^ a p (iu+ jv+kw) has vanishing 
divergence. Thus, for any steady flow, a vector potential 
iAi+JA 2 + KA 3 exists that is defined by the equation 
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or. 


Pi 

In Cartesian coordinates, 

Pu = 

Pv = 
pw = 


• curl A 

aA 3 

dA 2 

¥" 

■ sr 

dA x 

dA 3 

6z 

- 55" 

bA 2 


dx 

5y" 


(ID 


Therefore, for every three-dimensional steady flow the following 
functions exist: 


V s 


s 



Differentiation of these expressions yields equations (10), provided 
that 


5 r ^2 d p ^-2 . / 10 

55 J 5z" 4y “ Sz J S5~ dy (12 

The foregoing argument suffices to show the existence of the func- 
tions and cp for three-dimensional steady flows in general. 
Requirement (12) is certainly met hy most flows encountered in 
practice . 

Apparently, the theorem quoted atr the beginning of the pre- 
ceding paragraph may be extended to provide that the three-dimensional 
vector potential may be written as any one of the three possible two- 
component vector potentials. The existence of each of the three 
arrangements may be contingent on a restriction similar to equa- 
tion (12), which, it might be noted, has been shown here to be suf- 
ficient but has not b een shown to be a necessary condition. 
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Equation (3d) is satisfied in two-dimensional (axially symmetric) 
steady flow by defining f such that 


„ . at 

pru = 55F 


prv m 




(13) 


and in three-dimensional flow by defining t and 9 such that 


pru ; ^ 


prv 


&|( „ 1 d 9 
5x r 5s 


> (14) 



J 


For two-dimensional flow, definitions (14) reduce to definitions (13) . 
The existence of ,the functions t and 9 defined in equations (I4j 
my be established in a manner similar to that employed to show the 
existence of the analogous functions appearing in equations (10) . 

The quantities t and 9 are hereinafter referred to as "com- 
ponents of the vector potential" because (in the Cartesian system, 
for example) equations (10 ) my formally be obtained from equa- 
tions (11) by setting Aj = y , A-j_ = - cp , Ag = 0. 

Differential equations for vector potential. - 

(a) Cartesian coordinates: Substituting equations (10) into 

equations (2) and (4) for steady flow and adopting the subscript 
notation for partial differentiation yield 
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- (Vx^z) 




-Px(x,z) + 



(15a) 


% 




= - p z (x,z) + 



^y 0 x “ (^x +c Pz) 0 y + ‘Py^z 


(15b) 


” p ^y p x + p ^y p z + P a 


(p^)y + _(p CP y)y_ 


[ + ppt^y) 


(15c) 


or, for Pr = 1 and zero beat transfer, equation (15c) is replaced by 


c p T + 2 


"(?' *s) 2 + (f’x) 


= constant 


(15d) 


(b) Implicit coordinates: Substituting definitions (14) into 

equations (3) and an equation analogous to equation (4) yields, for 
steady flow. 
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r % 



x 


1 

r 


Vy 




EL A- Cp CD 

r 2 pr y y 


(16a) 


. El All cp 
+ r 2 pr ? y 


- - Pg(x,s) + 




?Vx -1(^4 ^) 


e y + 4 Vs 

r -<2 


pr ^yPx + pr 2 ^y®s + 4 


(1613) 


Pl_ * _( plr<Py )y 



+ A(pe y ) 
Pr y y 

(16c) 


or, for Pr = 1 and zero heat transfer, equation (16c) is replaced by 

constant (l6d) 


C P T + 2 


Ci ? v +(-^-9 y = 

ypr yy + \pr Jj _ 


Boundary conditions on vector potential. - The usual Boundary con- 
ditions on velocities applied in problems of the steady laminar bound- 
ary layer are, using Cartesian coordinates. 


■gjz.Q2 
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u(x Jtx »z) = u x (x,z) or t y (x,a,,z) 
w(x,<*,z) = w (x,z) or <? (x /CD ,z) 


PlUi(x,z) 

P^Cx^z) 


u(x,0,z) = 0 or t y.(x,D,z) = 0 
w(x,0,z) = 0 or 9 y (x,0,z) = 0 
v(x,0,z) = 0 or ty x (x,0,z) + <P z (x,0,z) * 0 


> (17) 


(18) 


An.oth.er boundary condition is usually required; for example, at the 
leading edge of an airfoil it is necessary to establish the condi- 
tions under which the initial growth of the boundary layer takes 
place. 


Boundary condition (18) , which involves a combination of ^ 
and cp , would probably be rather awkward to apply; the separation 
of V and with respect to their boundary conditions therefore 
seems desirable. The following argument shows that this separation 
can be made, so that condition (18) may be replaced by 

y (x,0,z) = cp(x,0,z) * 0 (19) 

Suppose that a solution has been obtained for a given problem, sub- 
ject to conditions (17) and (18). Equation (18) implies that 


y x (x,0,z) * hi(x,z) 
cp z (x,0,z) * h x (x,z) 

where hi is some arbitrary function. Thus, 


y (x,0,z) = J' hq dx + 1 x 2 ( 2 ) 

cp(x,0,z) = - I hi dz + h 3 (x) 


l * 

where h2 and h 3 are arbitrary functions. New quantities ¥ 
and cp* are now defined: 
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V* = ¥ - j h x dx - h ? (z) 

9 * = Cp+ J^ ^ dz " h 3 (x) 
Thus, in accordance with definitions (10), 


V*y =^y = pu 

- ~V*z = “*x -*z = PV 
9* y = cp y = pw 

The functions ty* and 9 * yield the correct velocity compo- 
nents, satisfy the same differential equation as ¥ and 9 , and sat- 
isfy conditions (17) and condition (19), which is a special case of 
condition (18). It is therefore correct to replace condition (18) 
ty condition (19), with the sole effect that the solutions for ¥ 
and cp will he made unique. 

The same argument applies in cases wherein definitions (14) 
are used.. 


Transformations of Equations of Motion 
Tiscosity -temperature relation. - It is assumed that the equation 


JL = JL 0 (20) 

Pod T qq 

may he used to represent adequately the variation of viscosity in the 
boundary layer. The constant C and the reference state T„> may he 
chosen to give the hest possible agreement with, for example, Suther- 
land’s formula, over the temperature range contemplated. A complete 
discussion of equation (20), as applied to flow with vanishing pres- 
sure gradient, is provided by Chapman and Rubesin in reference 10. 
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Howarth's transformations. - From reference 2, the following 
transformations are made in order to "bring equations (15) and (16) 
into forms approaching those of the corresponding incompressible 
equations : 



I (21) 


Before transforming equations (15) and (16) according to defini- 
tions (21) , it is desirable to make the physical quantities involved 
dimensionless. Hereinafter, the following quantities on the left 
will be considered to be measured relative to the quantities on the 
rigiht, where the subscript » denotes some reference condition: 


u,v,w relative to Ugo 

x»Y* z,r,X,Y,Z relative to UoP/p^u^ 


P relative to Poo 
T relative to v^/2c-g 

O 

p relative to p^u^ 

\jr, qi (as in equations (1©) and (15)) relative to 
V, cp (as in equations (14) and (16)) relative to 
0 relative to u ^ 


>( 22 ) 


(M 6D C/p 00 tloo)Q»««D 

(p ( _ c C/ft»Uo 3 ) 2 ft»u fc 
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With the conventions (22) and equation (20) , equation (15a), 
for example, "becomes 



(23) 


From equations (21), inasmuch as p is not to he considered a func- 
tion of y. 


i-ir” i. 

b _ b . bY b 

3x ~ 3x + Sx 5 y 

d d bY d 

V mn ” VMM —L. ■ l . l 

dz dZ. dz 5 y 

and from equations (10), 


u -*J 


W = ~CPy =9 Y J 

After relations (24) are introduced, equation (23) becomes 


(24) 


(25) 


(- 


p %( \x +Y x\y^ ' p ^% +c fe. +Y x% +Y z (: PY + 1 + | cpj 'hY + 


p Cpy ^1 Z +Y z%y) = ' PX + P Y (f 


ZfOZ 
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or, using the equation of state (2f), equation (26a) is obtained. The 
rest of equations (26), corresponding to equations (15) in Cartesian 
coordinates, and equations (28), corresponding to equations (16) in 
implicit coordinates, are obtained in a similar manner. 

In Cartesian coordinates, 

■*Y^nr - + - - i px + \ (j^*)**! (26a) 

+ + < 2Sb) 

" (^c + %) 0 y + < Py 0 z 

+ fF 0 YY (26o) 


iCfYPy-PyPz) + | 


(^Yy)^ + C5yy) ' 


JL , 

P» 


T_ 

Too 


(26d) 


or, for Pr = 1 and zero heat transfer, equation (26c) is replaced by 


T + O^y) 2 + ( c fy) 2 = constant 


(26e) 


Combining equations (17), (19j_, and (25) leads to the following boun- 
dary conditions on V and <P : 


Vy(X,C0,Z) * v^fXtZ) 
<P y(X/»,Z) « W^(X,Z) 


> (26f) 


ty(X,0,Z) =cp Y (X,0,Z) « Tjf (X,0,Z) =qp(X,0,Z) = 0 


Boundary conditions on temperature are required (only at Y * eo if 
equation (26e) is used) and a condition must be imposed to describe 
the initial growth of the boundary layers . (See Boundary conditions 
on vector potential. ) 
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In implicit coordinates, from equations (14), 


ru - - 7y 


rw » i q, y = q> y 


Equations (16) become 




r* ,m \?. 


Y s r Z' T' 


- ' l * + I I*) (F*) x 




Vi%j (»*J r 


1 (®Xs.l I 1 ST 


■-?%*! rt?’ 1?% +9 


| *x 8 x-(+ x 4 i ’s > 0 y‘F c( ’ a 


£?(^y%^;Ps) + 5F vf IT v®Y + ^ (%) 2+ ( 9 Yy) + SF 0 yy 


I = pi 

Poo T ~ (28d 

or, for Pr = 1 and zero heat transfer, equation (28c) is replaced by 


T + JL I (^•)^+( < Py) 2 I = constant 
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Two advantages are provided by use of transformation (21): 

1. The velocities are directly related to the vector potential; 
thus the satisfaction of boundary conditions is easier than under the 
original formulation, which included a variation with density. 

2. The transformed equations themselves closely resemble the 
corresponding equations of incompressible motion. 

Mangier 1 s transformations. - Equations (27) and (28) maybe 
brought into forms approaching those of equations (25) and (26) by 
means of a transformation introduced by Mangier (reference 8) for the 
purpose of relating axially symmetric flow to plane flow. The fol- 
lowing transformation is to be applied: 

'•X 

| s k 2 r 2 dX 

q z krY 
£ = s 
E s k ? 

G = k qT 

■where k is an arbitrary constant. The functions p, P, and T are 
not transformed. Thus, 



a 

3x 



a a 
5T + n x 



a 


d * i- 

3s 3? 
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U . i? Y - ^ 


v = - <?y = Gtj 


Under definitions (29), equation (28a) transforms as follows: 


rF^(k 2 r 2 F^+n x F T)Tj) ' (kr%£ )rkF nn + - r*(G n ) 2 


- 7*$>i * I ( A3 T *4*) 




which may be divided throughout "by k 2 r^, yielding equation (31a) . 
The rest of equations (31) are obtained in a similar manner: 


F r) F T) £ " F TH * k 2 r 3 S F n i " k2 r 3 


? p e + I l'p F+ ^3 t?* G j + F mn 


F n°n t ~ ( F e + ^ H) 


+ p^3 t + ^3 


^hhi + I ( t f+ ^ 1 ^ G ) % + 


(31a) 


(31b) 


G n e ^ 


i ( v 4 ^ 


2 Ip k 2 r 3 


F n 0 ^ - ^ 1+^3 e n + k 2 r 3 G n e i 
V ~i G ) e n + } ( FtjP | + ^3 Gt » p c) + (F nn } +(G n 


^ ) + e nn 


(31c) 
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JL = p 

Poo 



(Sid) 


or, for Pr 


1 and zero heat transfer, equation (31c) is replaced by 


T+(F^) 2 +(G^) 2 = constant 


(31e) 


Equations (27) and the analog of condition (19) lead to the following 
boundary conditions on F and H: 


Frj ( ) = n^(^C) 
Gt]( ~ w l( ^ ) 


> (31f) 


^(1,0,6) - G n (£,0,0 « F(e,0,0 = &(t,0,$) = 0 


Boundary conditions on temperature are needed (only at q = a> if 
equation (31e) is used) and the initial growth of the boundary layer 
must be described. (See Boundary conditions on vector potential . ) 

Transformation (29) thus confers the following advantages:, 

1. The velocities are related to the potential in the same 
manner as in Cartesian coordinates, thus simplifying the application 
of boundary conditions. 

2. The transformed equations approach Cartesian form. 

Further direct simplification of the boundary -layer equations 
does not seem feasible. Further simplifications, especially with 
regard to the reduction of the number of independent variables, must 
be sought in the consideration of special cases or classes of problems. 


APPLICATIONS 

Reduction to Problems of Known Solution 

The following examples are chosen to show how the foregoing 
theory specializes to some of the known boundary -layer flows. 
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Two-dimensional flow. - If it is supposed that two-dimensional 
flow obtains in the x,y -plane, then 

w = 0 


a 

a? 



whence , from equations ( 25 ) , 


u 

W - 0 


and, for example, equations (26) become 


__ i i Pt _ 

^Y^XY " VyY ' P P X + 2 T ^YY + \yY 


i y * + | ye v + (^)“ + 


*Y 0 X " V X°Y p T Y^X ' 2 J 


e. 


> (32) 


YY' f? YY 


and so forth. These are the equations of motion for plane flow, in 
terms of the stream function \j/ . 

Axially symmetric flow. - In the case of axially symmetric flow, 
the implicit coordinate system is used and r(x) is identified with 
the radius of a body of revolution. Thus, by the axial symmetry, 

w = 0 


a a 

5s 51 


= o 


whence, from equations (30), 


u = E 




G * 0 


and equations (31) reduce to equations of the same form as equa- 
tions (32) . 
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Flov about yawed infinite cylinder. - The Cartesian coordinate 
system, is used for flow about a yawed infinite cylinder; the x,y-plane 
is taken as the cross-sectional plane of the cylinder and the coor- 
dinate z is measured parallel to the axis of the cylinder. Thus, 
the potential flow has components in hoth the x- and z -directions. 
Therefore, in this case. 




0 


w £ 0 

qp / 0 


and equations (26) become 


1 ^X— — . - 


^Y^XY ‘ %% " P P X + 2 p t^YY + %YY 

W © - \1/ © = — ^¥m + m 

' Y TY f X y YY 2 V T£ y YYY 


L (33) 


% 9 x - Vi- }' *1% + I'S + & 9 yy 


If Pr = 1 and the wall temperature is constant, then, from equa- 
tions (8) and (26e), 

T + (¥y) 2 + (9y) 2 = a + bf Y 


where a and b are constants to be determined from boundary con- 
ditions on the temperature profile. 

Because of compressibility, the analytical separation of \jTy and 
9 y discussed in references 5 to 7 does not strictly apply. If incom- 
pressible flow is considered, then equations (33) become 
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-f PX+ + YYY 
V#XY '^X%Y = Cp YYY 

where p is constant. It is then possible to solve first for ^ , 
and subsequently to solve for , knowing % 


Boundary Layer on Flat Plates of Arbitrary Leading-Edge Contour 

The three-dimensional body considered here and the coordinates 
used in the analysis of 'its effect on a uniform stream are shown in 
figure 1(a). The surface of the body is a flat plate at zero angle 
of attack with respect to the stream. 


The differential equations of motion are, from equations (26a) 
and (26b), written for zero pressure gradient as 

Vrr - + 5y*zy ’ »yyy 

VfPxY ' + ^#ZY = ^YYY (34b) 


where 

u =?Y 
w = 

When u = ut = constant, the boundary conditions are 

^ y(X,os,Z) = ^y(Xq,Y,Z) — 1 
y(X,oo,Z) ~ *iPy(XQ,Y,Z) — 0 
Fy(X,0,Z) = \Rx,0,Z) » 0 
¥ y (x,o,z) =9(X,0,Z) = 0 


(34c) 


(35a) 
(35b) 
(35c ) 
(35d) 
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The result 9=0 obviously satisfies equation (34b) and boun- 
dary conditions (35b) and (35d) . This solution is presumably unique 
and may be interpreted to mean that in this case there is no secondary 
flow in the absence of a pressure gradient transverse to the stream. 
Equation (34a) then becomes 

Vy^XY " ^X^YY = ^YYY 


The solution of equation (36) , subject to boundary condi- 
tions (35a) and (35c), may be obtained by defining, by analogy with 
plane flow over a flat plate (reference 1, paragraph 53), 


Jx-Xq(z) 

and 

f = >v/X-Xq &(°) 


(37) 


(38a) 


whence , 


u = \j/ y = 6' (a) (38b) 

The application of equations (37 ) and (38) transforms equation (36) 
as follows : 


£ 0 ' ' + 20 ' ' * = 0 

The appropriate boundary conditions are 

0’ (on) = 1 

0* (O) = 0(0) - 0 


(39) 


(40) 


The solution of equation (59) subject to conditions (40) was given 
by Blasius and is presented in reference 1 (paragraph 53). 

Clearly, this solution to the problem of figure 1(a) is valid 
for any form Xo(Z) of the leading-edge contour. In particular, it 
is valid for the special configuration shown in figure 1(b). In this 
case, however, the Z-derivatives of flow properties in the boundary 
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layer (for example , which is related to the shear stress pu z ) 
suffer discontinuities at the plane Z = 0. This discontinuity may 
he shown by writing, using equation (37), 


do Xo'(z) 
3z 2 x-x 0 (z) 


5a 


and noting that Xq 1 (Z) discontinuously changes sign at the plane 
Z = 0. Not only is this circumstance physically inadmissable. 


but it violates the boundary -layer assumption that 


a 2 >:> a 2 


dY 2 dZ 2 


Apparently, there exists a narrow region extending in the stream 
direction from the apex of the body in figure 1(b) wherein the 
usual boundary -layer assumptions, and hence equations (33), do not 
apply. Outside of this region, hereinafter referred to as "wake", 
the usual quasi -two-dimensional solution of the boundary -layer 
equations for an infinite yawed flat plate is correct. 


A qualitative indication of the existence of such a wake in 
the flow shown in figure 1(b) is obtained by imagining a viscous 
wake in the form of a parabolic cylinder extending downstream from' 
each point of the leading edge. On either side of the X,Y-plane, 
these wakes then provide parabolic envelopes that, although not tan- 
gent to each other, are each tangent to the cylindrical wake pro- 
ceeding from the origin. Thus, in the neighborhood of the X,Y-plane, 
where the parabolic envelopes join to the cylindrical wake, it is 
clearly incorrect to make the usual boundary -layer assumption 



5Z 2 5y 2 


It is natural to suppose that the "wake" shaded region in fig- 
ure 1(b) has an effective width of the order of the boundary -layer 
thickness and thus that the gross effects of viscosity on the body 
are adequately given by the solution of equations (39) and (40). 
There is the possibility, of course, that the wake herein discussed 
has an over -all effect on the flow due to a difference in stability 
characteristics, as compared with the surrounding boundary layer. 


A wake such as is discussed in the preceding paragraph extends 
downstream from any corner in a leading-edge contour and presents a 
situation similar to that existing close to the leading edge of the 
flat plate, irrespective of the shape of the leading edge; as is 


V^iCYJ 
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well known, the "boundary -layer equations do not apply in the immediate 
vicinity of the leading edge. 

The inapplicability of the boundary -layer assumptions to a wake 
region extending downstream from a leading-edge comer is also asso • 
ciated with body configurations neighboring that of the flat plate 
with a leading-edge comer; that is, if the plate of figure 1(b) is 
given a small thickness, the solution of the boundary -layer equations 
then shows a variation of u^ at the plane Z = 0 which, though 

not discontinuous, is too rapid to be consistent with the boundary - 

3 2 d 2 

layer assumption — =■ » — = . As the thickness is further increased, 

dY 2 

this variation becomes less and less rapid, and, for some order of 
thickness, becomes consistent with the boundary -layer assumption. A 
similar argument applies if the curvature of the leading edge is 
imagined to be changed from infinity at the comer to some large 
finite value. 


Boundary Layer Associated with Supersonic Conical 
Potential Flows 

If the inviscid equations of motion are subject to boundary con- 
ditions on fluid properties that are constant along rays having a 
common focal point, it has been shown (reference 11) that the solution 
yields fluid properties that are constant along any such ray in the 
flow field; that is, that fluid properties (velocity, pressure, and 
so forth) are constant along each of a family of rays proceeding from 
a common focal point. This property is referred to as "conical 
symmetry." Solutions of this type exist, in general, only for super- 
sonic flows. 

For boundary -layer calculations, the significant feature of 
flows wherein the outer (potential) solution is conical is that pres- 
sure and velocity gradients at the outer edge of the boundary layer 
vanish in the direction of rays from a focal point. In this circum- 
stance, the equations of motion (31) appropriate to the implicit 
coordinate system are used. Distance along a ray is denoted by x, 
distance normal to the ray in the surface of the body by xs (that 
is, r(x) - x) 9 and distance normal to the surface by y. (See 
fig. 2(a).) Equations (29) thus become 
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rx 


l =k 2 


x 2 dX = |k Z X 3 


(41) 


11 * kxY - KXY 
£ = S 

Equations (31) may therefore "be written, inasmuch as, hy conical 
symmetry of the potential flow, p^ = 0, 

GF + F 
p tjtj nnn 


'n t ‘ ( F e + h a (j 'in + 5C Vnj - sj ( % )2 " Ff IL 

| ” ( F t + ^ G j) + sj a i°'i j * <Vi 


(42a) 


M 


3fe- p ’ < £ ) + ^ El p aG<5 'n +0 i 

- ( F e + h°i + * 


(421)) 


3lS 0 £ 


= 3^ S p '<0 + gf + < 


^qii) Z+ (^TiTi ^ 2 | T pp v t)T) 




e r 


(42c) 


The circumstance that p f = 0 suggests that the Blasius simil- 
“l/2 » 

arity variable T]£ (reference 1, paragraph 53) can he employed; 

that is, the boundary -layer development along rays from the origin 
may be expected to be parabolic. Thus, the following definitions 
are made, by analogy with the Blasius analysis: 
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X s n f 1 ^ 2 = 43 tl 1 ^ 2 
F = | 1/2 f (\>£ ) 

l/s , 

g “4 s(x>£) 

0 E 0(X,£) 
p = p(X,£) 

T = T(X,£) 

whence, 

u = Ftj = f x 

w = Gq = g^ 

Thus, equations (42) "become 


K (43) 


i f+i g+| g, ) f 




3 p °‘3 °£/ ~XX " 3 g X f X^ + 3 (g x } + 2f xxx 0 ^ 44a ^ 


/„lp'(t) 2 \ 2 

^ f+ 3 p s+ 3 S £/ S 'XX " 3 g X%£ 


I g x f X ■ h v ' { 0 + 2g xxx= 0 


(44h) 


/ i p’ ( d 2 \ 2 2 

\^ + 3 “p g+ 3 8.CJ S X - 3 g X e £ + 3P s X p ' + 


(f XX } +(g XX> 


+ 2 k - 0 


(44c) 


or, for Pr = 1 and zero heat transfer , equation (44c) is replaced "by 


T + (f X ) 2 + (g X ) 2 


= constant 


(44d) 
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or, using the equation of state, 

JL + (f, + (h = constant (44e) 

Poo P K A 

r 

c 

If 

The value of p given by equation (44e) may he substituted directly p 

into equations (44a) and (44h). The following boundary conditions 
apply : 

- ui( 0 

g v (»,£) = v ± U) 

f x (o,C) = g x (o,£) * f(o,0 * g(o, C) - o 

T(od, 0 = Tjl(£) (if equation (44d) is used) 

The equations (44) involve (if equation (44e) is used) two 
dependent variables f and g and two independent variables A. 
and £ . It seems reasonable to conclude that the solution of these 
equations for many cases is feasible provided modem high-speed com- 
puting techniques are used. 

In the light of the discussion of flat-plate flow contained in 
the preceding section, certain observations may be made concerning 
applications of equations (44) : 

The flow about the body shown in figure 1(b) may' be considered 
a special case of boundary -layer flow about a body with conical 
symmetry. In figure 2(b) this flat delta wing is drawn with reference 
to the coordinate system used herein for conical bodies. It may be 
shown that the solution given in the preceding section in Cartesian 
coordinates is consistent with the form required by equations (43) 
in conical coordinates: From equation (38b), the velocity com- 

ponents in conical coordinates (fig. 2(b)) are seen to be 

u = ( q) cos £ ~j 



w 


* - g* (o) sin £ 


( 45 ) 
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The distance X-Xq in Cartesian coordinates (fig* 1(b)) is, in con 
ical coordinates, 

X cos £ - X sin £ cot £o 

Thus, the a of equation (37) is, in conical coordinates. 


<7 * Y (X cos £ — X sin £ cot £ 0 ) 
- 1/2 


- 1/2 


YX 


cos C (1 - cot £ cot £ 0 ) 


- 1/2 


or, from equations (43), 


a = -i- k 

f5 


cos £ (1 - cot £ cot £ 0 ) 


- 1/2 


(46) 


Therefore, when equations (45) and (46) are combined, it is seen that 

u = u( \ £) 


w = w(\,£) 

as supposed in equations (43). 

The example just discussed concerns one of a class of conical 
bodies for which the boundary layer contains a wake of the type dis- 
cussed in the preceding section; that is, the boundary layer contains 
a narrow region to which the boundary -layer assumptions, and hence 
equations (43), do not apply. This class of bodies includes flat 
(or nearly flat) delta wings or rectangular wing tips. 


CONCLUSIONS 

Certain of the difficulties encountered in dealing with the equa- 
tions of three-dimensional motion for a laminar compressible boundary 
layer may be overcome by: 

1. The introduction of the vector potential in a manner that 
permits the expression of the three velocity components in terms of 
two scalar functions . 
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2. The use of a transformation of coordinates that has the 
simplifying effect of relating the vector potential to the velocity 
components in the same manner as for incompressible flow. Further, 
the transformed equations are of nearly incompressible form. This 
transformation requires the use of a linear type of viscosity- 
temperature relation. 

3. The use of a further transformation of coordinates in cases 
where an implicit coordinate system is employed (for example, when 
axial or conical symmetry is involved), such that the velocity com- 
ponents and vector potential are related in the same way as in Car- 
tesian coordinates. The transformed equations are of nearly Carte - 
sian form. 

The condition of constant total enthalpy is a solution to the 
energy equation in three- as well as in two-dimensional cases, when 
the Prandtl number is'l and for zero heat transfer. When the 
restriction of zero heat transfer is removed, an additional solution 
exists for the total enthalpy depending linearly on the velocity 
component in the direction of isobars of the flow, provided this 
component is constant in the outer flow. 

The equations herein developed for three-dimensional motion 
reduce directly to two-dimensional form for plane or axially sym- 
metric boundary conditions. 


In cases of flow over flat plates at zero angle of attack, when 
the leading edge is some arbitrary curve, the flow viewed in planes 
parallel to the stream and perpendicular to the plate is given by 
the plane Blasius solution. When the leading edge has a corner, 
this solution contains discontinuous derivatives at the normal plane 
passing through the comer in the streamwise direction. It is inferred 
that , in a narrow region extending downstream from the comer , a 
"smoothing" process occurs to which the boundary -layer equations are 
not applicable. This effect is not expected to be physically impor- 
tant, although the possibility exists that it affects stability and 
transition. 

When the boundary layer is associated with supersonic potential 
flow having conical symmetry, the Blasius similarity variable may be 
applied with respect to variations in normal planes containing the 
apex. Thus, the laminar boundary -layer development is parabolic 
along rays from the apex. The use of this information reduces the 
number of independent variables to two and is thus considered to bring 



NACA TR 2279 


35 


such problems within the range of effectiveness of modem computing 
techniques, especially if the assumptions of Prandtl number equal 
to 1 and zero heat transfer are made. 


Lewis Flight Propulsion Laboratory, 

Rational Advisory Committee for Aeronautics, 
Cleveland, Ohio, September 5, 1950. 
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